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We consider a number operator-annihilation operator uncertainty as a well behaved alternative to 
the number-phase uncertainty relation, and examine its properties. We find a formulation in which 
the bound on the product of uncertainties depends on the expectation value of the particle number. 
Thus, while the bound is not a constant, it is a quantity that can easily be controlled in many 
systems. The uncertainty relation is approximately saturated by number-phase intelligent states. 
This allows us to define amplitude squeezing, connecting coherent states to Fock states, without a 
reference to a phase operator. We propose several setups for an experimental verification. 
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I. INTRODUCTION 

Finding a phase operator conjugate to the number op- 
erator and constructing number-phase uncertainty rela- 
tions has an extensive literature 1-3]. However, defining 
a Hermitian phase operator for an infinite system, i.e., a 
harmonic oscillator, is not possible, and all the different 
approaches must make certain compromises. 

Historically, the first important contribution was that 
of Dirac, who introduced the phase observable tj> based 
on the decomposition of the annihilation operator as a = 
Rexp(i<p) Q. Assuming that R and <j) are Hermitian 
operators, one obtains R — N 1 / 2 , where N = a^a is the 
number operator. Hence, exp(i</>) must be equal to 



E 



\n) (n+l| = (7V+l)- 1/2 a, 



(1) 



However, E is not unitary, thus <f> cannot be Hermitian 
either. 

Several methods have been presented to circumvent the 
difficulties above. Since there are extensive reviews on 
the topic [1H3J , we cite only the literature that is direct!; 
connected to our approach. Susskind and Glogower 
constructed the Hermitain operators C = ^(E + E^) and 
S = ^{E — 2?t) to describe the quantum phase. They 
obtained uncertainty relations with them, however, for 
the description of the phase two operators were needed. 
In order to overcome this inconvenience, Levy-Leblond 
Q suggested the use of the non-Hermitian E defined in 
Eq. ([T]). He argued that physical quantities could also 
be represented by non-Hermitian operators, interpreted 
the meaning of variance for such operators and wrote 
down uncertainty relations with N and E. Later, Her- 



mitianphase operators were constructed for finite sys- 
tems [3|, lZ|- This makes it possible to carry out a cal- 
culation for an expression with the phase operator for 
a finite dimension D, and then take the limit D — > 00, 
which provides the value corresponding to the infinite 
dimensional case. By use of this theoretical background, 
the number-phase uncertainty relation could be obtained 
and the states saturating them, called number-phase in- 
telligent states, were identified [8[- The procedure that 
provides a connection from one number-phase intelligent 
state to another one with a smaller number variance is 
called amplitude squeezing 

Connected to these ideas, in this paper we choose the 
two operators to be not N and E — (N + l) _1 / 2 a, but 



simply N and the annihilation operator a. 
the relation 



(AN) 2 



(Ao) 



> 



(N) 



We present 



(2) 



We will show that states saturating the number-phase 
uncertainty are very close to saturating Eq. ([2]). This 
makes it possible to define amplitude squeezing without 
a reference to a phase operator. 

In addition to its connections to quantum optics, this 
problem is also interesting from the point of view of quan- 
tum information theory. A family of uncertainty rela- 
tions with N and a has already appeared in Ref. pij] , and 
has been used for the detection of quantum entanglement 
pTH13| . Such uncertainty relations made it possible to 
construct entanglement conditions with small experimen- 
tal requirements. Remarkably, these conditions detect 
non-Gaussian entangled states that cannot be detected 
based on the first and second moments of the quadra- 
ture components TJ|. The uncertainty relation Eq. ^fy 
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presented in this paper can be seen as a single relation 
replacing the family of uncertainty relations described in 
Ref. For given (N), Eq. ([2} identifies most of the 

values for the variances of iV and a that are not allowed 
by quantum physics [l5j . 



2 



The paper is organized as follows. In Sec. II, we dis- 
cuss how the variance of the annihilation operator can 
be defined. In Sec. Ill, we derive the uncertainty relation 
Eq. @. In Sec. IV, we discuss the tightness of the uncer- 
tainty relation presented. Finally, in Sec. V we discuss 
possible physical tests of the proposed uncertainty rela- 
tion. In the Appendix we present uncertainty relations 
for two-mode systems. 



II. VARIANCE OF THE ANNIHILATION 
OPERATOR 

In this section, we will discuss the definition and the 
properties of the variance of the annihilation operator. 
We will relate it to quadrature independent properties of 
the quantum state. 

We define the variance of a non-Hermitian operator A 
as dEl] 



(AA) 2 = (A^A) - (A*) (A). 



(3) 



Note that, for non-Hermitian operators usually we have 
(AA) 2 ^ (AA^) 2 . 

Let us now consider the A = a case. (Aa) 2 is zero only 
for coherent states. The variance (Aa) 2 measures, in a 
sense, how close the quantum state is to a coherent state. 
For this reason, it has been used to study the dynamics 
of various quantum systems (e.g., see Refs. [l7H20l |). 

Let us now interpret (Aa) 2 by relating it to the quadra- 
ture components. 

(i) Let us define the quadrature components as 



xp 

P0 



ae +i P - at e - tf 



where f3 is real. Then, one finds that 
(Aa) 



2 _ (A^) 2 + (A^) 2 1 



(4) 



(5) 



Hence, (Axg) 2 + (App) 2 is independent of the choice of 
the angle /? 

When discussing the invariance properties of (Aa) 2 , it 
is instructive to point out its connection to the correlation 
matrix defined as 



■P = 

(Axp) 2 \{(AxpApp + AppAxp)) 

\{{Ax p App + AppAxp)) (Ap/j) 2 



Thus, since (Axp) 2 + (App) 2 is independent of (3, it 
seems to be a good measure of the uncertainty of the or- 
thogonal quadrature components. Note that an alterna- 
tive measure could be the product (Axp) 2 (App) 2 , how- 
ever, it is not independent from j3. 

(ii) In another context, (Aa) 2 can be expressed as 



Z7T 



2tt 



0=0 



[Axpfaip - - 



(7) 



Thus, (Aa) 2 is connected to the average variance of the 
quadrature components xp. That is, if /3 is chosen ran- 
domly between and 2ir according to a uniform proba- 
bility distribution then (Aa) 2 + | gives the expectation 
value of the quadrature variance (Aa:^) 2 . 

(iii) Finally, let us examine the connection between 
(Aa) 2 and important properties of the Wigner function of 
the quantum state. For the following discussion, as well 
as in the rest of the paper, we will leave the /3 subscript, 
and will use x and p in the sense of xq and po , respectively. 
(Aa) 2 gives information on the sharpness of the peak of 
the Wigner function W(x,p) of the state since [22| 



(Ax) 2 + (Apf 



[(x - (x) f + {p- (p)) 2 ] W(x,p)dxdp. 

(8) 



For states with a non- negative Wigner function (i.e., 
squeezed coherent states), 2(Aa) 2 + 1 is the sum of the 
squared widths of the Wigner function in two orthogonal 
directions. The sharpest peak is obtained for the coher- 
ent states for which Eq. (|8]) is the smallest. 



III. UNCERTAINTY RELATION WITH THE 
NUMBER AND THE ANNIHILATION 
OPERATORS 

In this section, we will present a simple derivation of 
Eq. <j2j) and relate it to the uncertainty relation with the 
variances of N and E. We will also discuss how to improve 
the relation Eq. J2|). 

We start from the two Heisenberg uncertainty relations 



(AiV) 2 (Ap) 2 > 
(AA0 2 (Aa;) 2 > 



\\{*)\\ 
\\{P)W 



(9) 



(6) 



where we used the fact that for operators A and B we 
have (AA) 2 (AB) 2 > \\([A,B])\ 2 Summing the two 
inequalities of Eq. ©, using Eq. ((SJ) and \ (p}\ 2 + \ (x)\ 2 — 
2|(a)| 2 , one obtains the following uncertainty relation 
with N and a 



One can obtain V pi from T p through orthogonal trans- 
formations. However, the trace of T, which equals to 
(Axp) 2 + (App) 2 , remains invariant under such transfor- 
mations. 



(AN) 2 



(Aa) 



>\\(a)\ 2 - 



(10) 



From Eq. (fTC))) it follows that knowing (a), which deter- 
mines the "center" of the Wigner function W(x,p) of 
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the state, and (Aa) 2 , which is based on the the width 
of the Wigner function in two orthogonal directions, a 
lower bound for the particle number fluctuation can be 
obtained. 

The bound in the uncertainty relation Eq. (|10[) is not 
a constant: it depends on (a), which is zero for a wide 
class of states. It would be meaningful to find a similar 
relation with a constant bound or at least with a bound 
depending on a quantity that is easily measurable and 
controllable. 

We will now construct a relation in which the bound 
depends on (N) rather than on (a). For that, we add 
[(Aa) 2 + ±]/4 to both sides of Eq. ((TDJ}, and using (Aa) 2 

Y - | (a) j 2 , we obtain Eq. ©. The right hand side of 
Eq. ((2J is minimal for the vacuum |0). In all other cases 
the right hand side is greater than | , thus the uncertainty 
finds some part of the (Aa) 2 -(AA^) 2 plane inaccessible 
for quantum states. 

Next, we will relate Eq. ([2} to the uncertainty relation 
with N and E. For that, we determine the form of Eq. © 
for the case of large N and (A7V) 2 <C (A^) 2 - In this case 



a s» y/ (N)E and we obtain 

(AA0 2 (A£) 2 > \ [1- (AS) 2 ] , 
which is in accordance with the results of Ref. 0] 



(11) 
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FIG. 1: (Color online) Numerical test of the inequality J2|. 
F refers to the Fock state [n = 25) , C refers to the coherent 
state \a = 5) , Z refers to the point that saturates inequality 
© for (Aa) 2 = 0. (solid) Boundary of the region defined by 
Eq. ([2]) for (TV) = 25. All points below this line correspond to 
aphysical (Aa) 2 ~(AN) 2 values, (dashed) Points correspond- 
ing to squeezed coherent states. The equation of the curve 
is given in Eq. (|23[) . Inset: (solid) Boundary of the region 
defined by Eq. (J2J. (dashed) Points corresponding to states 
with a Gaussian wave vector Eq. (I16|l . (circles) States corre- 
sponding to photon added coherent states Eq. $27} . 



(AN) 2 (AE) 2 > - |l - (AE) 2 - (P<°>) 



(12) 



where = |0) (0| . 

Finally, Eq. (|2|) can be improved by use of the 
Robertson- Schrodinger inequalities (22[. First, they can 
be used to improve the two uncertainty relations in 
Eq. (|9]). Then, after steps similar to the previous ones, 
we obtain 



(AA0 2 + 7 



(Aa) 



1 



-^ + rt l({AiV ' Aa}+>|2 ' 

(13) 



where {A, B} + = AB + BA is the anticommutator. 



IV. TIGHTNESS OF THE INEQUALITY 

In this section, we will investigate the tightness of 
Eq. ©, and will look for quantum states that are close 
to saturating it. We will also discuss that the states sat- 
urating the left hand side of Eq. © interpolate between 
coherent states and Fock states. 

Our inequality does not contain the highest possible 
lower bound. The reason for that is that we constructed 
Eq. ([2|) by summing the two uncertainty relations in 
Eq. (jQj) . While the new relation Eq. (fTO]) is valid, it is 
not tight since the two uncertainty relations in Eq. © 
are saturated by different states. Thus, the tightness of 
the bound in Eq. ^ must be verified. 

In Fig. 1, wc plotted the points corresponding to values 
of [(Aa) 2 , (AN) 2 ] that saturate Eq. © for (AT) = 25. All 



points below this line violate the relation Eq. ©. For 
Fock states, 

(Aa) 2 ock = (A0, (AAf)§ ock = 0. (14) 

Hence, Fock states saturate Eq. ©. For coherent states 
we have 



(Aa) c 2 oh = 0, 



(AiV) 2 oh 



(AO. 



(15) 



For (Aa) 2 = 0, the particle number variance saturat- 
ing Eq. © is (AN) 2 = h{N). This already shows that 
the lower bound in Eq. ([2f cannot be optimal because for 
(Aa) 2 = there is no quantum state with a smaller parti- 
cle number variance than (N). States minimizing (AA^) 2 
for < (Aa) 2 < (N), in a sense, interpolate between 
coherent states and Fock states. 

We will now examine the tightness of Eq. ([2]) numer- 
ically through choice of appropriate trial states. Our 
search can be simplified by noting that it is sufficient 
to search over wave vectors with non-negative real ele- 
ments. To see this, let us consider a state of the form 
\^>) = ^2 n \c n \e l ^ n \n) . One finds that if all angles <p n are 
set to zero, (AAf) 2 and (AT) do not change. On the other 
hand, \(a} \ cannot decrease. Hence, (Aa) 2 = (N) — |(a)| 2 
cannot increase. Thus, it is sufficient to search over states 
with all 4> n = 0. Moreover, a state with <fi n = const, x n 
will give the same values for (AN) 2 ,(N) and (Aa) 2 as 
does a state with d> n = 0. 
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FIG. 2: (Color online) Numerical test of the inequality 
Eq. ([2|. The difference between the left-hand side and the 
right hand-side of Eq. <(2j) divided by the right hand side is 
shown for states Eq. (|16[) with a Gaussian wave vector for 
(from left to right) No = 100, 25, 10 and 5. 




FIG. 3: (Color online) Distance of the points corresponding 
to states with a Gaussian wave vector Eq. (|16|) from the curve 
corresponding to states that saturate Eq. ((2| for (from left to 
right) iVo = 100, 25, 10 and 5. A and No are the parameters 
of the state Eq. (JTSJl . 



Gaussian wave vector 



Let us consider states with a Gaussian state vector 



(n - N f 



4A 2 



(16) 



where C is for normalization. For such states, (N) ~ Nq 
and (AN) 2 A 2 . In the inset of Fig. 1, the dashed 
line corresponds to states of the form Eq. (fTB"|) while the 
solid line corresponds to points saturating Eq. (j2]). It 
can be seen that states of the form Eq. (fT()|) are close 
to saturating Eq. @. In Fig. 2, we plotted the relative 
difference between the left hand side and the right hand 
side of Eq. for particular values of A and Nq. 

Finally, Fig. 3 shows the distance in the (Aa) 2 -(A7V) 2 
plane of the points corresponding to states of the form 
Eq. ([16]) from the curve corresponding to states that sat- 
urate Eq. 23]. The distance does not grow with N 
and remains smaller than 0.15. This means that for large 
N, for which (Aa) 2 and (AiV) 2 cannot be measured with 
an accuracy of 0.15, the states Eq. [[To]) are practically 
the intelligent states of the uncertainty relation Eq. ((5J). 

The states Eq. (fTo]) are the subset of the number-phase 
intelligent states called \g') presented in Ref. |8|. There, 
the coefficients of |n) have a Gaussian dependence on 
n, just as in Eq. (|16|) . however, the phase of the coeffi- 
cients is not zero but has a linear dependence on n. As 
we have already discussed, a state vector with a phase 
with a linear dependence on n has the same values for 
(Aa) 2 , (AJV) 2 and (N) as a state with zero phase. Thus, 
all the \g') states presented in Ref. [1] are very close to 
saturating Eq. ([2]). Hence, our inequality makes it possi- 
ble to define number-phase intelligent states and ampli- 
tude squeezing |9| without a reference to a phase opera- 
tor. 



There are other states known to be number-phase in- 
telligent state s l24l I25] ] . We now consider the states pre- 
sented in Ref. [26| . They are defined as the superposition 
of coherent states on a circle 



-f CO 



\ao,u) oc / exp(— ^u 2 4> 2 — i8<j)) |cvoe^) c 
J — 00 

where 8 = a 2 ,. The overlap with Fock states is 
a™ 

(ao,u\n) oc — p= exp 



(17) 



(n-5f 



= (oio\n) x exp 



2u 2 
(n~6) 2 
2u 2 



(18) 



The second expression stresses the fact that we have the 
overlap of a coherent state \cto) and a Fock state |n) , 
multiplied by a Gaussian centered around a 2 ,, that is, the 
expectation value of the particle number for |cvo) . Thus, 
\ao, u) has an almost Gaussian wave vector for large N in 
the number basis. Hence, these states give similar results 
numerically for our uncertainty relation Eq. ([2]) as does 
Eq. CEU). 



B. Squeezed coherent states 

It is natural to ask to what extent squeezed coherent 
states approach the curve defined by Eq. ©. Squeezed 
coherent states can be obtained from the vacuum state 
as MM 



\a,0=D(a)S{Q\0), 



(19) 



where D is the displacement operator and S is the squeez- 
ing operator. Next, we will use the relationships 

D^(a)aD(a) = a + a, 
S\C)aS{() = acosh(s)-a t e i,? sinh(s), (20) 
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where £ — se 



i i) 



(A«)L 



l«,0 

2 



Hence, with a = \a\e , we obtain 
sinh 2 (s) + \a\ 2 , 
sinh 2 s, 



(AN) 2 aX) = H 2 [cosh(2s)-sinh(2s)cos(26>-tf)] 



2 sinh 2 (s) [l +sinh 2 (s)] 



(21) 



For given |a| and s, the variance (AiV) 2 in Eq. ([21]) is 
minimal if cos(2# — d) = 1. This is fulfilled, for example, if 
9 = i? = 0, that is, both £ and a are real and nonnegative. 
Hence, 

(AN) 2 min (\a\,s) = H 2 [cosh(2 S ) - sinh(2 S )] 

+ 2 sinh 2 (s) [1 + sinh 2 (s)] . (22) 

Based on Eq. (|22|). we obtain the smallest possible (AiV) 2 
for squeezed coherent states, for given (Aa) 2 and (N) as 

(AAOLn = [(N) - (Aa) 2 ] + (Aa) 2 - ^Aaj\ 



2(Aa) 2 [l + (Aa) 2 



(23) 



A dashed curve in Fig. 1 corresponds to Eq. 

Let us interpret this result. Since (Aa) 2 and (AiV) 2 
are invariant under a rotation around the origin in the 
x-p plane, we can start from coherent states |a) with a 
real and positive a. Then, the state we considered for the 
curve Eq. (|2"3")l corresponds to squeezing of the x quadra- 
ture component, which is called "number-squeezing" in 
the literature (e.g., see Ref. [28j]) and it reduces the num- 
ber variance for a small amount of squeezing. Thus, for 
small squeezing Eq. (|23p is not far from the bound given 
by Eq. ([2]). With further squeezing, the number vari- 
ance starts to grow. Thus, for large (Aa) 2 there are 
no squeezed coherent states giving an almost minimal 
particle number variance, and one has to look for non- 
Gaussian states for that. As a by-product of our discus- 
sion, note that the non-Gaussianness of quantum states 
can be verified by measuring only (Aa) 2 and (AiV) 2 . 



C. Displaced Fock states 

Displaced Fock states are defined as (29| 
\a, n) — D(a) \n) . 
Using Eq. (|2U|) . we obtain 

( N )\a,n) = n + \a\ 2 , 
(AN) 2 a . n) = (2n+l)|a| 2 , 
(Aa)f Qin) = n. (25) 
Hence, for displaced Fock states we get the equation 

(AiV) 2 = [2(Aa) 2 + 1] [(N) - (Aa) 2 ] , (26) 

where (Aa) 2 must be a non-negative integer. It is fulfilled 
by both Fock states and coherent states. Other points in 
the (Aa) 2 — (AiV) 2 plane satisfying Eq. $Zjo§ are very far 
from saturating Eq. @. 



D. Photon-added coherent states 

Photon added coherent states are defined as [3(| 

\a,m) oc (a)) m \a) . (27) 

They are close to saturating Eq. (|2|), as can be seen in 
Fig. 1. 

E. Eigenstates of a^a + const, x a 

According to Heisenberg's method, states that mini- 
mize the uncertainty product (AX) 2 (AY) 2 for Hermitian 
X and Y with a constant commutator are the eigenstates 
of X+icY, where c is some constant [3l| . While in Eq. ^ 
we do not have the product of the uncertainties of two 
Hermitian observables, this method can still give us the 
idea of considering the states \d,k) defined through the 
eigenvalue equation 



(a) a + da) \d, k) = k \d, k) 



(28) 



where d and k are constants. Writing \d, k) as J2k c 'n \ n ) i 
we obtain 



(k 



--n+l 



dy/n+1 



(29) 



for the coefficients. Equation ([29]) leads to a normalizable 
wave vector only if k is a non-negative integer. In this 
case, q = for all I > k. Numerical evidence suggests 
that states |d, A;) are close to saturating Eq. ©, but they 
are inferior to the states given by Eq. (fTS)) . 



F. States minimizing (AiV) 2 for given (Aa) 2 and (N) 

Let us now look for states that minimize (AN) 2 for 
given (Aa) 2 and (N). For that, we will follow an approach 
similar to the one presented in Ref. 31|. Since (Aa) 2 = 
(N) - \((x) 2 + (p) 2 ), this task can be reformulated as a 



search for the states that minimize (AiV) for given 
(p), and (N). Let us write the state as |$) = c n 



(24) 

Hence, we have to look for the minimum of the function 

/(($|,|<i>),A J v,A p ,A ;c ) = (iV 2 )-iV 2 

+X x ((x) - x ) + X p ((p) -po) + X N ((N)-N ), 

(30) 

where are Lagrange multipliers. Note that we do not 
include explicitly the ( < I ) | < I ) ) = 1 condition in the function 
/. The minimum is given by one of the critical points 
for which all derivatives are zero. Equation Q30p can be 
rewritten as 



/«$|,|$),Ajv,Ap,A x ) = (0(X N ,X p ,X x )) m , 



(31) 
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where O is defined as 



0(X N , X p , X x ) 



Ajva^a + {a) a) 2 + 



V2 



X x — iAj; 

~7T 



(32) 



We have to look for , A x fe) , X p k) ) that mini- 

mize Eq. (|3"Tj) . It is easy to see that |$( fc )^ must mini- 
mize (0(X ( ^,\ ( p\\ x k) )}. Hence, states must be 

the eigenstates of the operator 0(X^\ X^ k \ Xx) with 
the smallest eigenvalue (i.e., they have to be "ground 
states"). Note that the operator given in Eq. (1521 . 
appears as a system Hamiltonian in self-consistent cal- 
culations for the Bose-Hubbard model based on the 
Gutzwiller ansatz [H,|33j]. 



V. DISCUSSION 

Let us discuss the use of quantum states that mini- 
mize (AiV) 2 for given (N) and (Aa) 2 . They present a 
trade-off between two requirements: the smallest possi- 
ble variance of a randomly chosen quadrature component 
and the smallest possible particle number variance. In a 
sense, they are similar to states minimizing (AN) 2 for- 
given (A(p) 2 . The latter present a trade-off between the 
smallest possible variances for phase measurements and 
for particle number measurements. 

Clearly, the right hand side of Eq. is not a constant, 
but it is a quantity that can be controlled easily in many 
systems. Moreover, note that the measurement of (Aa) 2 
does not require the measurement of variances of x and 
p if we use (Aa) 2 = (N) - |(a)| 2 = (N) - ±((x) 2 + (p) 2 ). 

A single trapped ion seems to be a good candidate for 
testing our inequalities and realizing quantum states that 
saturate them [34T - l36j ] . For a trapped ion, x and p are 
the physical position and momentum coordinates, and N 
determines the energy of the ion. 

The uncertainty relation Eq. ^ can also be verified ex- 
perimentally in a single mode electromagnetic field. The 
two orthogonal quadrature components can be measured 
for example with homodyne detection [33] ■ The result is 
not influenced by which two orthogonal components we 
choose to measure. 

Bose-Einstein condensates of alkali atoms seem to be 
also a possible candidate for experiments [3^, 38]. It is 
usual to talk about number squeezing in multi-well Bose- 
Einstein condensates in the sense that an increase in the 
the barrier hight between the wells decreases the number 
fluctuation within the wells j39|. Here, one has to note 
that for cold atoms the particle number is conserved. Be- 
cause of that, for a single bosonic mode of cold atoms, su- 
perpositions of states with different particle numbers are 
not allowed. For this reason, for pure states (AN) 2 = 0. 
It is possible to mix states with different particle num- 
bers making (AN) 2 > 0. However, even for such states 



(a) = and (Aa) 2 = (N), which makes our inequalities 
trivial in such systems. 

While we cannot create particles in a single mode, 
we can move particles from one mode to another one. 
Thus, it is instructive to consider two-mode systems of 
cold atoms. The two modes can be realized with atoms 
in a double well or with a single Bose-Einstein conden- 
sate of two state atoms. Let us denote the annihilation 
operators of the two modes by a\ and a2, respectively. 
The corresponding particle numbers are N\ and N%. If 
(Ni) < (N 2 ) and (AN 2 ) 2 < (N 2 ) 2 then with the substi- 
tution 



a — > 



N 



Ni, 



(33) 

the uncertainty relations Eqs. (flOl) and ([2]) can be tested. 
In the Appendix we present relations that do not require 
such approximations. 

Finally, in the statistical physics of bosonic systems, 
(fy(x)}, i.e., the expectation value of the field operator 
plays the role of the order parameter. In this context, 
our findings present a quantitative relationship between 
the variance of the field operator and the variance of the 
particle density ^(x)^(x). 



VI. SUMMARY 

We constructed uncertainty relations with the particle 
number and the annihilation operator. The variance of 
the latter describes the uncertainty in the phase space, 
and is independent of the absolute phase of the quadra- 
ture components. We proposed quantum optical systems 
in which our inequality could be tested. 
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Appendix: Uncertainty relations for two-mode 

systems 

For the two-mode system, inequalities similar to 
Eqs. (fTU)) and ^ can be found using the Schwinger rep- 
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resentations of the angular momentum operators 

T 



1 



01 



a i 

«2 



(A.l) 



for I = x,y, z where 07 are the Pauli spin matrices. Let 
us define an operator that is an analog of a in the two 
mode-system as 



(A.2) 



a = J x — ij y = a\a 2 - 

With this definition, we have 

|AS| 2 ) = i[(AS) 2 + (ASt) 2 ]^(AJ ;c ) 2 + (AJ 2/ ) 2 , 

(A.3) 



2 _ / t \2 , / t \2 

y) ■ 



\{a)V = (J*Y + (J, 



Using Eq. (IA.3I) and the Heisenberg uncertainty rela- 
tion (AJfc) 2 (AJ;) 2 > j|(Jm)| 2 ! we obtain the analog of 
Eq. (TO} 



(A7V 1 ) 2 (|A5| 2 )> i |(5)| 2 . 



(A.4) 



Adding ±(|Aa| 2 ) to both sides of Eq. (|A.4j) and using 



Jl + Jy = ~(Ni + l)(N 2 + l)--, (A.5) 



we obtain an analogue of Eq. ^ 



(AAq) 5 



|Aa| 2 )>-((iY 1 + l)(7V 2 + l))--. 



(A.6) 



As we mentioned previously, number squeezing with 
Bose-Einstein condensates in a double-well can occur if 
the barrier between the wells increases [39] . Equation. ([2]) 
bounds the number variance of a well in such systems 

Si. 
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